clc;
clear;
R = 220;
C = 33e-12;
L = 100e-6;

omega0 = 1/sqrt(L*C);

% Резонансная частота
f0 = omega0 / 2 / pi;
fprintf('f0 = %f MHz\n', f0/1e6); 

% Ось времени
T = 1 / (10000*f0);
t = 0:T:(15 * 1/f0);
lt = length(t); 


% Внешние параметры, внешние фазовые переменные
std=3;


E=std*randn(1,length(t));


 

% Начальное состояние
Uc(1) = 0;
% Начальное напряжение конденсатора
i(1) = 0;
% В начальный момент ток в цепи отсутствует
di(1) = 0;
Ul(1) = di(1) * L;
Ur(1) = di(1) * R;
% Напряжение на катушке

for k = 2:lt

Ul(k) = E(k) - Uc(k-1) - i(k-1)*R;
% Напряжение на катушке
di(k) = Ul(k) / L;
% вызывает приращение тока,
i(k) = i(k-1) + di(k-1)*T;
dUc(k) = i(k) / C;
% который определяет заряд конденсатора
Uc(k) = Uc(k-1) + dUc(k)*T;
% Напряжение на резисторе
Ur(k) = i(k-1)*R + di(k-1)*T*R;
U = Ur;
end 

figure(1);
plot(t*1e9, [Ur; E; U]);
xlabel('t, ns')
ylabel('U_r, E, U, Volt')
legend('U(t)', 'E(t)', 'U(t)')
grid on

std_U = sqrt(cov(U));
fprintf('std_U = %std_U B'); 

fmin=0e12;
fmax=5e6;
A=1;
f=fmin:((fmax-fmin)/100):fmax;
%K=nan(1,length(f));
i=sqrt(-1);  
    k = (i.*f*2*pi*R*C)./(-1.*f.*f*4*pi^2*L*C+1+i.*f*2*pi*R*C);
    for l=1:length(f)
    b(l)=abs(k(l));
    end
    
    
figure(2);
plot(f, b);




